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Abstract 

We consider heat conduction across an ordered oscillator chain with harmonic 
interparticle interactions and also onsite harmonic potentials. The onsite spring con- 
stant is the same for all sites excepting the boundary sites. The chain is connected 
to Ohmic heat reservoirs at different temperatures. We use an approach following 
from a direct solution of the Langevin equations of motion. This works both in the 
classical and quantum regimes. In the classical case we obtain an exact formula for 
the heat current in the limit of system size N — > oo. In special cases this reduces 
to earlier results obtained by Rieder, Lebowitz and Lieb and by Nakazawa. We also 
obtain results for the quantum mechanical case where we study the temperature 
dependence of the heat current. We briefly discuss results in higher dimensions. 
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1 Introduction 



Rieder, Lebowitz and Lieb [T] (RLL) considered the problem of heat conduction across 
a one- dimensional ordered harmonic chain connected to stochastic heat baths at the two 
ends. The main results of this paper were: (i) the temperature in the bulk of the sys- 
tem was a constant equal to the mean of the two bath temperatures, (ii) the heat cur- 
rent approaches a constant value for large system sizes and an exact expression for this 
was obtained. RLL considered the case where only interparticle potentials were present. 
Nakazawa [2] (N) extended these results to the case with a constant onsite harmonic 
potential at all sites and also to higher dimensions. 

The approach followed in both the RLL and N papers was to obtain the exact nonequi- 
librium stationary state measure which, for this quadratic problem, is a Gaussian distri- 
bution. A complete solution for the correlation matrix was obtained and from this one 
could obtain both the steady state temperature profile and the heat current. 

In this paper we use a different formalism to calculate the heat current in ordered 
harmonic lattices connected to Ohmic reservoirs (for a classical system this is white noise 
Langevin dynamics). The formalism has been discussed in detail in [31 S] and follows from 
a direct solution, by Fourier transforms, of the Langevin equations of motion. A general 
formal expression for the heat current can be obtained and this has the form of an integral 
of the heat transmitted over all frequencies. For the one-dimensional case this expression 
for current was first obtained by Casher and Lebowitz [3] using different methods. An 
advantage of the approach used here is that it can be easily generalized to the quantum 
mechanical regime [U EJ [7J [8j [9]. As shown in [I] the results obtained are identical to 
those obtained using the nonequilibrium Green's function method [TUl [TT] . Here we show 
how exact expressions for the asymptotic current (N — > oo) can be obtained from this 
approach. We also briefly discuss the model in the quantum regime and extensions to 
higher dimensions. 

The model we consider here is slightly different from the Nakazawa model. We consider 
the pinning potentials at the boundary sites to be different from the bulk sites. This allows 
us to obtain both the RLL and N results as limiting cases. Also it seems that this model 
more closely mimics the experimental situation. In experiments the boundary sites would 
be interacting with fixed reservoirs which can be modeled by an effective spring constant 
that is expected to be different from the interparticle spring constant in the bulk. We also 
note here that the constant onsite potential present along the wire relates to experimental 
situations such as that of heat transport in a molecular wire attached to a substrate or, 
in the two-dimensional case, a monolayer on a substrate. Another example would be the 
heat current contribution from the optical modes of a polar crystal. 

The paper is organized as follows. In SecfS] we introduce the model and derive the 
main results of the paper. In Secj3] we present results for the quantum mechanical case 
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and in SecJH generalizations to the higher dimensional case. Finally we conclude with a 
short discussion in Sec El 

2 Model and results in the classical case 




Figure 1: A schematic description of the model. 

We consider N particles of equal masses m connected to each other by harmonic springs 
of equal spring constants k. The particles are also pinned by onsite quadratic potentials 
with strengths k Q at all sites except the boundary sites where the pinning strengths are 
k Q + k! [see Fig. (pQ)]. The Hamiltonian of the harmonic chain is thus: 

N 1 1 N ~ X 1 1 

H = Z^™^ + 2 k ° X ^ + ^ 2^ Xl+1 ~ Xl)2 + 2 k '^ 2 + ^ ' (2,1) 
i=i i=i 

where xi denotes the displacement of the particle at site I from its equilibrium position. 
The particles 1 and iV at the two ends are immersed in heat baths at temperature Tl 
and Tr respectively. The heat baths are assumed to be modeled by Langevin equations 
corresponding to Ohmic baths. In the classical case the steady state heat current from 
left to right reservoir is given by [31 E] : 

Jc = MTL ~ TR) f du T N (u), (2.2) 

where Tn(u) = 7 2 u; 2 |GiAr| 2 , 

G = \-muj 2 I + $ -S]" 1 , 

®im = {k + k' + k ) 8i }m - k 8i tTn -i for 1 = 1, 

= -k 5i >m+1 + (2k + k Q ) 5i, m - k <5 Zjm _i for 2 < I < N - 1 , 

= (k + k' + k ) 8i tm - k 8i, m+1 for I = N , 

T.im = i^u8i m [8n + 8i N ] , 
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and J is a unit matrix. We now write G = Z~ l /k, where Z is a tri-diagonal matrix 
with Zu = Z NN = {k + k Q + k' — muj 2 — ijuj)/k, all other diagonal elements equal 
to 2 + k D /k — muj 2 /k and all off-diagonal elements equal to — f . Then it can be shown 
easily that |Giat(w)| = l/(k |Ajv|) where An is the determinant of the matrix Z. This is 
straightforward to obtain and after some rearrangements we get: 

A at = [a(g) sin Nq + 6(g) cos Nq]/ sing , (2.3) 

7 2 cj 2 k' 2 2k\ 2k' 2i-fu t ,k! 

where a(q) = [2 --— + — -—] cosg + — - 2 - -—[1 + (-- 1) cosg] , 

rv rv rv rv rv rv 

n 2 u 2 k' 2 2k' . 2iju.iV ^ . 

6 (<?) = [-^- TJ + ^r} sm q + ^-(j-i)wq, 

f\ IX IX fx rv 

and q is given by the relation 2k cos q = —muo 2 + k + 2k. This relation implies that for 
frequencies outside the phonon band k < muo 2 < k Q + 2k the wavevector q becomes imag- 
inary and hence the transmission coefficient T(u) decays exponentially with N. Hence 
for large N we need only consider the range < q < ix and the current is given by: 



2j 2 k B (T L -T R ) r^duj 



UO 



Jc = ' T,_ ' / dq\^\ ^ , (2.4) 



k 2 7T J dq \A 
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with muo 2 = k Q + 2k[l — cos (g)]. Now we state the following result: 

where gi(q) and (72(9) are an Y t wo well-behaved functions. This result can be proved 
by making an expansion of the factor 1/[1 + g 2 (q) sin (Nq)] (valid for \g\ < 1 in the 
integration range), taking the N —>■ 00 limit and resumming the resulting series. Noting 
now that An can be written as |Aat| 2 = (|a| 2 + |6| 2 )[1 + r sin(2Nq + </>)]/ [2 sin 2 (g)] where 
rcos0 = (ab* + a*b)/(\a\ 2 + \b\ 2 ), rsin0 = (|6| 2 - \a\ 2 )/(\a\ 2 + |6| 2 ), we see that Eq. f[2T4"j) 
has the same structure as the left hand side of Eq. (12.51) . Hence using Eq. (12.51) and after 
some simplification, we finally get: 



Ji 



c 



jk 2 k B (T L - T R ) r sin 2 qdq 



nm J A — ficosg 

yk 2 k B (T L - T R 



mQ 2 



(A - VA 2 - n 2 ) , (2.6) 



where A = 2k(k - k') + k' z + '— and Q = 2k(k - k') + 



m m 
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Eq. (12. 6p is the central result of this paper. We now show that two different special 
cases lead to the RLL and N results. First in the case of fixed ends and without onsite 
potentials, i.e. k' — k and k Q = 0, we recover the RLL result []]: 



jRll _ kk B (T L -T R ) 



c ~ 27 



v v / 4 

1 + \ 1 + - 

2 2 V v 



mk 

where v = — — . (2.7) 
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The can also be obtained using the RLL approach [12] and agrees with 

the result in Eq. (12. 61) . In the other case of free ends, i.e. k' = 0, we get the N result [2]: 



N _ k>yk B (T L - T R ) 
c ~ 2(mfc + 7 2 ) 



, A A 4 
1 + - - 



where A = - ^ 9N . (2.8) 



2 2V A J k{mk + 7 2 



3 Quantum mechanical case 

In the quantum case the heat current across a chain described by the Hamiltonian Eq. (12. ip 
and connected to Ohmic heat baths is given by [I]: 

1 [°° 

J Q = - du hwT N (u)[f(uj,T L ) - f(u,T R )] , (3.1) 

where f(u,T) = l/[e^/( fcsT ) — 1] is the phonon distribution function and Tn is as given 
in Eq. (12.21) . Here we consider the linear response regime where the applied temperature 
difference AT = Tl — Tr << T with T = (Tl+Tr) /2. Expanding the phonon distribution 
functions f(u),TL 7 R) about the mean temperature T we get the following expression for 
the current: 

Jc .MZ^£^(-_y cosech2( -_ )rM . ( , 2) 

We then proceed through the same asymptotic analysis as in the previous section and get, 
in the limit iV — ► oo: 

-yk 2 h 2 (T L -T R ) r , sm 2 q 2 2 ( fvuj q 



^ = 7 L l R) / ^ A "I" ^ 2 cosech 2 f^), (3.3) 

Q A-KknmT 2 L y A-ficoso q V2fc R T/' v ; 



AnkBrnT 2 J A — Vlcosq 
where u 2 = [k Q + 2k(l — cos q)]/m . 

We are not able to perform the above integral exactly. Numerically it is easy to obtain 
the integral for given parameter values and here we examine the temperature dependence 
of the current (note that in the classical case the current depends only on the temperature 
difference). In Fig. (2) we plot the current as a function of temperature in three different 
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cases (i) k' = k,k Q = 0, (ii) k' = 0, k Q = and (iii) k' = 0, k Q ^ 0. Particularly interesting 
is the low temperature (T << hik/rn) 1 ' 2 /ks) behaviour (shown in inset of Fig.(2)) which 
is very different for the three cases. The low temperature behaviour can be obtained 
analytically by examining the integrand at small q. We then find for the three different 
cases: (i) J Q ~ T 3 , (ii) J Q ~ T and (iii) J Q ~ e -^°/( fcflT ) /T 1 / 2 , where w = (fc^m) 1 / 2 . 




Figure 2: Plot of the scaled heat current with temperature (in units of fiikjwb) 1 ! 2 /ks) for 
three different parameter regimes (see text). Inset shows the low temperature behaviour. 



4 Higher dimensions 

Heat conduction in ordered harmonic lattices in more than one dimension was first con- 
sidered by Nakazawa [2]. The problem can be reduced to an effectively one-dimensional 
problem. For the sake of completeness we reproduce their arguments here and also give 
the quantum generalization. 
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Let us consider a d- dimensional hypercubic lattice with lattice sites labelled by the 
vector 1 = {l a },a = l,2...d, where each l a takes values from 1 to L a . The total number 
of lattice sites is thus N = L\L2...Ld- We assume that heat conduction takes place in 
the a — d direction. Periodic boundary conditions are imposed in the remaining d — 1 
transverse directions. The Hamiltonian is described by a scalar displacement X\ and as in 
the ID case we consider nearest neighbour harmonic interactions with a spring constant 
k and harmonic onsite pinning at all sites with spring constant k Q . All boundary particles 
at Id = 1 and Id = Ld are additionally pinned by harmonic springs with stiffness k' and 
follow Langevin dynamics corresponding to baths at temperatures Tl and Tr respectively. 

Let us write 1 = (lt,ld) where L, = (h,l2---ld-i)- Also let q = (qi,q2---Qd-i) with 
q a = 2im/L a where n goes from 1 to L a . Then defining variables 

X M = T y 2 i rl/2 E , (4-1) 

one finds that, for each fixed q, Xi d (q) {l d = l,2...L d ) satisfy Langevin equations corre- 
sponding to the ID Hamiltonian in Eq. (12. ip with the onsite spring constant k Q replaced 
by 

A(q) = k a + 2[d - 1 - cos (^)l- ( 4 - 2 ) 

a=l,d-l 

For Ld — > oo, the heat current J(q) for each mode with given q is then simply given by 
Eq. (l2.6p with k Q replaced by A q . In the quantum mechanical case we use Eq. (13.31) . The 
heat current per bond is then given by: 

J = LL \ E J (q)- ( 4 - 3 ) 

Note that the result holds for finite lengths in the transverse direction. For infinite trans- 
verse lengths we get J = f ... f Q n dqJ(q) / (2n) d ^ 1 . 



5 Summary 

In this paper we have derived the exact formula for the heat current through an ordered 
harmonic chain in the limit of infinite system size. Our derivation is different from the 
methods used by RLL [1] and N [2] and is for a slightly different version of the models 
studied by them. We have presented the quantum mechanical generalization of the results. 
In that case one gets, in the linear response regime, a temperature dependent current with 
interesting low-temperature behaviour. We have also stated the results for the general 
(^-dimensional case. 
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